Magnetoelectric coupling in a ferroelectric/ferromagnetic chain revealed by 

ferromagnetic resonance 
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Understanding the multiferroic coupling is one of the key issues in the field of multiferroics. As 
shown here theoretically, the ferromagnetic resonance (FMR) renders possible an access to the 
magnetoelectric coupling coefficient in composite multiferroics. This we evidence by a detailed 
analysis and numerical calculations of FMR in an unstrained chain of BaTiOs in the tetragonal 
phase in contact with Fe, including the effect of depolarizing field. The spectra of the absorbed 
power in FMR are found to be sensitive to the orientation of the interface electric polarization and 
to an applied static electric field. Here we propose a method for measuring the magnetoelectric 
coupling coefficient by means of FMR. 
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Introduction.- Materials with multi ferroic (magnetic, 
electric, and/or elastic) orders, called multiferroics (MP), 
have attracted increased attention again [11-01, niainly 
due to the discovery that the notoriously small multifer- 
roic coupling in bulk matter may well be increased by 
a controlled engineering of low dimensional compounds, 
opening thus the way for the design of qualitatively new 
device concepts 0-01 • For instance, magnetoelectricity 
allows for the control of magnetism with an electric field 
which has a large potential for environmentally friendly 
sensorics and spintronics applications with low-energy 
consumption. The progress in this field depends criti- 
cally on the understanding of the magnetoelectric (ME) 
coupling and on developing methods to assess its prop- 
erties. This is particularly so, as currently several ME 
coupl ing mechanisms arc being discussed, e.g. in Refs. 
Isl-fToj . On the other hand, an established approach 
for probing the ferromagnetic response is the ferromag- 
netic resonance in which the sample is usually 
subjected to crossed static and time-dependent magnetic 
fields. Hence, it is natural and timely to envisage a pos- 
sible mapping of the multiferroic dynamics in an FMR 
setup with the aim to draw conclusions on the nature 
of the ME coupling and relate it to the PMR signal. In 
fact, PMR has been experimentally shown [H, to be 
sensitive to acoustic waves in ferromagnetic-ferroelectric 
structure. To our knowledge, multiferroic PMR, as sug- 
gested below has not yet been realized experimentally for 
the chosen interface, though several studies are known for 
multiferroic interfaces with other types of ME-coupling 
[13, 13 • To be specific, we focus on a special class of 
ME coupled materials, so-called composite MPs [gI [19|. 
that may be synthesized from a wide range of mate- 
rials that, when composed together yield a strong ME 
coupling and stable ferroelectric (EE) and ferromagnetic 
(FM) orders at room temperatures. An example that 
has been studied intensively, theoretically and experi- 
mentally is BaTi03(BT0)/Pe The ME coupHng 



in this system is predicted to be an interfacial effect and 
relatively high l2ll. [28l . Isil l32j . An experimental evidence 
is presented in [20]. 

For a reliable prediction, a realistic modelling is indis- 
pensable since the ME coupling is relatively weak in to- 
tal (due to its interfacial nature) and the MF dynamics 
is governed by a series of interrelated effects, as shown 
below. The route followed here is based on a combina- 
tion of the Landau-Lifshits-Gilbert and the Ginzburg- 
Landau dynamics using the total MF free energy Fs 
density [13, HI, US, including the ferroelectric Ffe, the 
ferromagnetic i^FM , and the part Sint that involves the 
ME interface coupling. FpE includes the energy densi- 
ties corresponding to the Ginzburg-Landau-Devonshire 
potential Fqld [la, [s^ , to the spatial inhomogcneity of 
the order parameter Fqe [37j , the depolarizing field con- 
tribution Fdf[31|, the dipole-dipole interaction i^^Dj, as 
well as the applied external electric field Faef- ^'fm 
[sot incorporates the nearest-neighbor exchange interac- 
tion i^EXC, the anisotropy contribution Fani, the FM 
dipole-dipole interaction F^^^ and the energy density of 
the interaction with an external magnetic field _Famf- 
Realistic parameters corresponding to bulk BaTiOa in 
the tetragonal phase and bulk Fe are employed. Our 
simulation show that the FMR spectra of the absorbed 
power are indeed sensitive to the FE order as controlled 
by an applied static electric field, evidencing thus the ac- 
cess of FMR to the ME coupling. The key finding of 
this study is a practical proposal to estimate the mag- 
netoelectric coupling coefficient from the FMR spectra. 
With the refinement in the spatiotemporal resolution of 
the FMR technique we expect the multiferroic FMR to 
play a vital role in uncovering the details of ME coupling. 
Theoretical model.- Our treatment is based on the 
Ginzburg-Landau phenomenology, i.e. we study the dy- 
namics of coarse-grained order parameters (P, and Mj ) of 
the FE/FM chain, that result from an averaging of the 
microscopic quantities over an appropriate cell. These 
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cells form in our calculations the sites i or j for the local 
Pi or Mj . At the interface (site i = j = 1) of the FE/FM 
composite the mobile spin-polarized electrons in the FM 
rearr ang e as to screen the electric polarization at the FE 
part [401 , leading thus in effect to a local ME coupling of 
the interface magnetization Mi with the interface polar- 
ization Pi that can be expressed [l^ as i?iNT = ^Pi ■ , 
where A is the ME coupling coefficient. The FE polariza- 
tion vector Pi develops in time according to the Landau- 
Khalatnikov (LKh) equation f4l| 



dP^ 
dt 



5F^ 

' SP, 



(1) 



with the FE relaxation constant = 2.5 • 

10~^ [Vms/C] [42|. The magnetization dynamics Mj 
is governed by the Landau-Lifshitz-Gilbert (LLG) 
equationjH, li^ 
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where 7 is the gyromagnetic ratio and apM is a Gilbert 
damping constant. The local effective field is defined as 
H]-^^{t) = ~f§^- The FMR power absorbed by the 
chain we infer from (cf. e.g., (45l. l46j) 



dt 



dt, (3) 



where fio is the vacuum magnetic constant, ofm is the cell 
size, A'^T and T = 27r/a; are the number and the period 
of the external magnetic field cycles, respectively. The 
sum runs over magnetization sites j. The total magnetic 
field applied to the system is H^{t) = He^ + Ho coswicy, 
whereas Hq <C H. 

It the following the imaginary part of the transverse 
magnetic susceptibility y" will be calculated, for which 
x"-^I^l^holds[ll,[ll. 

Contact of ultrathin FE and FM (E = Q).- First, we nu- 
merically model a contact of a single FE site and a single 
FM site in zero E-ficld. Employing the FE free energy 
density for the tetragonal phase of BaTiOsjlJ] 
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and the FM free energy density for a uniaxial crystal in 
an external magnetic field [l^l 
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(5) 



we calculated the spectra shown in Fig. f. The posi- 
tion of the peak for the case of zero ME coupling follows 
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FIG. 1. FMR spectra for the multiferroic chain of A''fe ~ 1 
and A''fm = 1-sites. BaTiOa has the tetragonal phase with 
the following coefficients in the free energy density (eq. (HV) 
apE = -2.77 • 10'' [Vm/Cljsl, /3fe = 1.7 ■ 10* [VmVC^][33 
and Ps = 0.265 [C/m^][33. The parameters of bulk Fe 
assumed in the calculations (eq. ([5])) for this figure are: 
Ki = 4.8 • lO-* [J/m^][3i, Ms = 1.71 • 10*^ [A/m][3il and 
QfFM = 0.1. For the time-dependent magnetic field it is taken 
HoHo = 28 • 10"=* [T], uj/{2tv) = 4 ■ lO'-* [Hz]. Inset shows the 
dependence of resonance peaks on A. 



according to which we find yUo^^rcs(A = 0) « 0.087 [T] 
(Fig. 1) [48[. For a finite ME coupling constant the reso- 
nance field is modified by the anisotropy induced by the 
screened polarization 
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This relation reveals that the resonance condition de- 
pends on the orientation of the surface polarization. In 
the simulations the spectra are calculated at equilibrium 
with the effect of opposite orientation of the polarization 
and magnetization. For a finite A the magnetization is 
aligned along z and the polarization Pz is negative. This 
results in a shift of the resonance peaks towards smaller 
fields (Fig. 1, A 7^ 0). The shift of resonance fields re- 
mains linear as a function of A (inset of Fig. 1). 
We also note, that the ME coupling does not affect the 
intensity and the width of the spectra, similarly to the 
effect of the uniaxial anisotropy axis oriented along z axis 
in macro-spin FM nanoparticles(46j . 
Contact of ultrathin FE and FM (E ^ Q).- Considering 
the MF chain in the presence of an electric field E which 
acts directly on the electrically active part of the chain, 
wc expect for i5 || z that the resonances shift to weaker 
magnetic fields (eq. (jG))). To favor the stability of the 
magnetization along the z direction, we increase the fre- 
quency of the oscillating magnetic field which lowers the 
intensity of x" ■ Similarly to Fig. 1, a finite A leads to a 
shift of the resonance field towards smaller fields (A = 
is not shown in Fig. 2), this is also predicted by a very 
recent study [2^. The nonzero electric field applied par- 
allel to the MF chain acts on the FE polarization and 
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FIG. 2. FMR spectra for the multiferroic chain of A''fe ~ 1 
and A'^PM = 1-sites. All parameters are adopted from Fig. 
1, except the ME coupling parameter A = 0.2 [s/F] and the 
frequency of the applied mag netic field w/(27r) = 6-10" [Hz]. 

indirectly changes the ME coupling energy enlarging the 
resonance fields (eq. ^ for the field E ~ Ee^). 
Let us define the shift of the resonance position from eq. 
([5]) for a finite electric field with a reference to zero E-field 

MO AH = if res ^ 0) - if res = 0) = 

AAP, = A[Pz(^^0)-Pz(S = 0)], (7) 

and plot Ho AH against APz (Fig. 3, left), we can deter- 
mine from the slope of the obtained dependence the ME 
coupling parameter A = 0.202 [s/F], the original value 
of which was set to 0.2 [s/F]. The analysis of the po- 
larization states shows that the abrupt increase of the 
resonance field for i? = 5 ■ 10^ [V/m] (Fig. 2) refers to 
the reversal of the polarization from the antiparallel to 
the parallel arrangement with respect to the magnetiza- 
tion orientation in the FM part. Since re-polarization 
of the FE modifies significantly the system and creates 
an energy jump at the interface responsible for a parallel 
orientation of FE and FM sites, it is natural to expect 
that the corresponding resonance point will deviate from 
the linear dependence, which is proved by Fig. 3 (left 
panel). Therefore, a linear fit of the spectral data should 
be performed for the points obtained for the electric fields 
that do not result in re-polarization of the ferroelectrics 
(E < 5 ■ 10^ [V/m]). If the spectral points measured for 
high fields were included into the fitting data set, the 
resulting values of A can exceed the actual value con- 
siderably. Notwithstanding the difficulties of FE polar- 
ization measurements [HI during the FMR experiment, 
the proposed procedure gives a transparent method for 
obtaining the magnitude of the ME-coupling coeSicient 

Thin FE/FM contact for E ^ 0.- For a thicker MF 
system the total FE energy density additionally includes 
interactions of the neighboring sites [l^ and the dipolc- 

dipole interactions Pj^giH pN>i = (^^f^P^ + 
^P4 + ^(P,,+i - P,,)2 - P,,£;,) + . The FM 



FIG. 3. Shifts of resonance positions calculated according to 
expression (O for different applied electric fields. Values of 
FE polarizations are calculated for the resonance magnetic 
fields presented in Fig. 2. Resonance curves for the fields 
E = {1.0; 3.0; 4.0; 4.5} • lO'' [V/m] are not shown in Fig. 2. 



energy density (eq. ([5])) is supplemented by the ex- 
change interaction and the FM dipole-dipole interaction 

^^^^iHi F^>^ = (Mz.f - 

fioMj ■ H^{t)j + P™i. The ME coupling in our case is 

limited to the vicinity of the FE/FM interface. Hence, 
we expect the infiucncc of ME coupling on MF dynamics 
to be less pronounced as compared with the single-state 
case (Fig. 2), which is evidenced by Fig. 4. Never- 
theless, the spectral lines are clearly distinguished for 
FE switching from antiparallel to parallel orientation re- 
garding the direction of magnetization vectors in FM 
part. The FE switching occurs for lower field values 
in comparison with single-state case because the pres- 
ence of interaction between the sites lowers the barrier 
of the Ginzburg-Landau-Devonshire potential. Applying 
the procedure outlined above, we plotted (inset to Fig. 
4) the variation of peak position versus the averaged fer- 
roelectric polarization. Remarkably, the obtained points 
fits well to a straight line for the case of low polariza- 
tion values corresponding to a non-switched ferroelectric 
layer. The points related to the switched value diverge 
from a linear scaling and were not considered. The slope 
of the linear fit obtained from the plot is 0.01376, which 
is lower than the value of A = 0.2 [s/F]. This effect most 
probably appears because the magnetoelectric coupling is 
limited to a single interface and then propagates through 
a chain of five sites at each side of the device, resulting in 
a smaller variation of FMR spectra. A straightforward 
step consisting in a multiplication of the fitted value of 
A by the number of sites yields the value of 0.0688 [s/F] 
that is three times smaller than A. We think that it is 
possible to find a proper re-normalization constant from 
geometric considerations, which, however, go beyond the 
scope of this letter. The most important result is that the 
peak position / polarization variation plots present the 
same linear dependence as that observed for the single- 
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FIG. 4. FMR spectra for a chain of A*'fe — 5 and A''fm ~ 5 
sites, cell size: apE = cifm = 5 ■ 10~® [m]. The applied field 
frequency is uj/{2it) = 30 • 10^ [Hz], FM exchange stiffness is 
A = 2.1-10"" [J/m] [si], and A = 0.2 [s/F]. Other parameters 
are as in Fig. 1. The FE coupling strength is calculated as 
KFE = Gii/alE, where Gn = 51. • 10"" [JmVC^]!!!. Insets 
show the peak on an enlarged scale resolving the states with 
non-switched and switched ferroelectric layer. The linear fit 
of the resonance shifts involves the non-switched polarization. 



site system, which endorse the proposed method for the 
measurement of the magnetoelectric coupling for thicker 
composite multiferroic systems. 

Remarks and conclusions.- As shown above, the FMR 
spectra of the absorbed power of an unstrained thin com- 
posite mutliferroic chain depend critically on the ME cou- 
pling. The peaks of resonance absorption are sensitive to 
the magnitude and the orientation of the FE polariza- 
tion vector (Fig. 1) in the absence of an electric field. 
Applying a static electric field changes the value of 



and causes a shift in the peak position according to (eq. 
([6])). Figs. 2 and 3 demonstrate how in our case the ME 
coupling can be accessed. At first, the static electric field 
is applied along the direction of the FE minima such that 
it shifts sizably the peaks relative to the field-free case. 
The shifts of the fields fioAH (eq. ^ are plotted against 
the measured FE polarization, yielding the ME coupling 
coefficient as the slope of the /xoAiJ(APz) dependence 
(Fig. 3, left panel). The plot oi plqAH{E) (Fig. 3, right 
panel) is nonlinear, hence, the ME coupling can only be 
determined for known P^iE) dependence. We note, the 
IJ.oAH{E) dependence will be different for the rhombo- 
hedral phase|33 of the FE 

Based on the results in (Figs. 3, 4), on the estimates 
made for the ME coupling coefficient 12J, l34i . and on the 
calculations for the elongated MF chains [33 we suggest 
to choose a MF contact consisting of a thin FM layer 
(electrode) and a thick FE. The thin FM layer serves to 
avoid additional peaks in the FMR spectra due to spin 
waves as well as to avoid a decay of the ME coupling 
in long FMs, whereas the thick FE part stabilizes the 
FE polarization for measuring P^ (Fig. 3, left). Note, for 
highlighting the ME dynamics dominated by the interfa- 
cial screening charges, the principal axis of Fe is rotated 
by 45° with respected to BTO[100], we have then a good 
lattice match (~ 1.4%) for the epitaxial growth of bcc Fe 
on the tetragonal BTO. 
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